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ABSTRACT 

The  main  aim  of  this  work  is  to  study  electron  scattering  in  imperfect  semiconductor 
heterostructurcs.  The  source  of  unintentional  disorder  is  the  interface  roughness  at  the 
hctcrojunctions  occurring  during  growth.  In  order  to  achieve  this  goal  W'e  solve  numeri- 
cally the  tw’o-dimensional  Ben  Danicl-Duke  equation  for  the  electron  scattering  problem. 
Our  model  assumes  open  boundary  conditions  along  the  growth  direction  and  periodic 
ones  parallel  to  the  heterojunctions.  We  then  compute  the  reflection  and  transmission  ma- 
trices that  govern  channel  mixing  due  to  interface  roughness  scattering.  The  knowdedge  of 
the  mixing  matrices  allow  us  to  calculate  the  transmission  coefficient  in  any  heterostruc- 
ture made  of  wide  gap  semiconductors.  As  an  example,  we  compute  the  transmission 
coefficient  in  resonant  tunneling  devices  based  on  double-barrier  structures. 

INTRODUCTION 

Electron  scattering  by  imperfect  heterojunctions  reduces  electron  mobility  due  to 
rough  surfaces  even  in  good-quality  heterostructures  [1].  So  far,  there  are  analytical 
results  concerning  the  propagation  of  wave  packets  in  a randomly  layered  medium  when 
the  potential  is  a random  function  of  only  one  coordinate  [2],  but  for  a small  number  of 
layers,  as  in  double  barrier  heterostructurcs  (DBH),  in-planc  disorder  becomes  important 
and  one  expects  such  approaches  to  fail.  Realistic  models  of  in-plane  disorder  usually  lead 
to  intractable  analytical  models;  hence  the  importance  of  numerically  solvable  models  to 
bridge  this  gap.  An  important  contribution  was  already  provided  by  Henrickson  et  al., 
who  applied  the  tight-binding  Green  function  method  to  account  for  electron  scattering 
by  interface  roughness  in  imperfect  DBH  [3].  How'ever,  in  this  paper  electron  transmission 
through  a DBH  w^as  described  by  a rather  artificial  model  of  disorder,  namely  periodic 
roughness  with  random  relative  phases  at  the  interfaces. 

In  the  present  w'ork  we  introduce  a two-dimensional  effective-mass  model  to  study  the 
effects  of  interface  roughness  scattering  on  electron  transmission  through  semiconductor 
heterostructures.  To  this  end,  the  Ben  Daniel-Duke  equation  is  discretized,  boundary 
conditions  are  discussed  and  scattering  solutions  are  found  by  means  of  the  transfer- 
matrix  method  for  any  arbitrary  heterostructure  made  of  wdde  gap  semiconductors.  The 
model  is  worked  out  in  a two-dimensional  space  for  computational  limitations,  although 
it  w'ill  be  clear  that  generalization  to  three  dimensions  is  rather  straightforw'ard.  Finally, 
W'e  present  the  numerical  results  for  the  transmission  coefficient  through  imperfect  DBH 
and  the  main  conclusions  of  the  work. 
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THEORY 


We  consider  the  Ben  Daniel-Duke  etiuation  with  constant  effective  mass  m*  at  the 
r valley  for  the  electron  envelope  function  z),  where  the  two  sjratial  directions  are 

shown  in  Fig.  1.  The  whole  structure  is  divided  into  three  different  regions,  namely  left 
(I)  and  right  (111)  contacts  and  the  DBH  (II),  where  scattering  by  lateral  disorder  takes 
place. 
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Figure  1.  Scln'inatic  view  of  the  sample.  Regions  I and  III  are  the  (dcu’tri- 
cal  contacts  of  the  sample  and  electrons  undergo  scattering  i)roeesses  by  lateral 
disorder  only  at  region  II  (DBH). 


We  consider  a mesh  with  lattice  spacings  Uy  and  o.-  in  the  y and  z directions,  re- 
spectively. Defining  ty  = -ti^ /(2m'u^)  and  L = we  obtain  the  discretized 

equation  for  the  envelope  function 

fi(0n  + l,m  d"  V’n  — l.m)  "h  t^(t'n.m4-l  "f  — l)  d'  (fbj.rn  2f~  (1) 

The  potential  term  U„,„,  in  Eq.  (1)  is  given  by  the  conduction-band  edge  energy  at 
the  i>oint  {nay.mad)  which,  in  turn,  depends  on  the  ,41  mole  fraction  in  the  vicinity  of 
that  position.  Therefore,  lateral  disorder  enters  the  (>(|uation  through  this  diagonal  tc'rm. 
Contacts  an'  characti'rized  by  fiat  band  conditions,  U{n,m)  = 0.  Solutions  of  Ecp  (1)  can 
be  determined  from  the  appropriate  boundary  conditions  by  using  tlu'  transh'r  matrix 
method.  The  boundary  conditions  are  open  in  the  z direction,  and  pi'riodic  on  each  slidt', 
that  is  in  the  y direction.  The  former  imply  plane  wave  .solutions  in  the  s axis,  and  the 
latter  yield  an  energy  discretization  on  y.  .4s  a consequence,  this  discretization  results 
in  a number  of  transverse  channels  equal  to  the  number  of  points  in  the  transverse  mesh 
direction.  Once  the  solutions  are  known,  the  transmission  coefficient  is  computed  as  a 
function  of  the  electron  energy. 
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MODEL  OF  DISORDER 


In  order  to  treat  lateral  disorder,  we  have  considered  the  occurrence  of  islands  on 
the  interface  between  two  consecutive  layers  having  identical  lateral  sizes  all  of  them  and 
being  consecutive  one  to  each  other.  In  our  model,  islands  have  heights  (measured  from 
mean  position  of  the  heterojunction)  that  are  randomly  distributed  to  mimic  random 
fluctuations  of  local  flux  of  atoms  during  growth.  It  is  then  feasible  to  express  the  rough 
profile  of  the  interface  between  two  consecutive  epilayers  defining  the  following  height 
function 

Kv)  = v^Wn{Q{y  - O + + i)C  -y]~  i}-  (2) 

n 

Here  h{y)  represents  the  deviation  from  the  flat  surface  at  position  y,  0 is  the  Heavyside 
theta  function,  C is  the  island  width,  w„  is  a random  variable  associated  to  the  n th  island 
that  controls  the  fluctuation  around  the  mean  value,  and  77  is  the  largest  deviation  — in 
absolute  value  --  assuming  that  the  w„’s  are  uniformly  distributed  between  —1  and  1. 
Hereafter  r)  will  be  referred  to  as  degree  of  lateral  disorder.  The  random  variables  Wn  take 
values  from  —1  to  1 uniformly,  satisfying  the  following  correlator  {wnWm}  — where 

the  brackets  indicate  average  over  different  realizations  of  the  disorder.  In  particular, 
notice  that  (h{rj))  = 0. 


Figure  2.  Islands  modeling  lateral  disorder  (roughness)  at  the  interface 
between  two  epilayers  (GaAs  and  Al^Gaa-xAS). 


RESULTS 

We  have  performed  several  numerical  calculations  in  order  to  study  the  effect  of  both 
lateral  and  compositional  disorder  over  the  transport  properties  of  DBH  made  of  GaAs- 
AlxGai_xAs  heterostructures.  We  start  by  considering  the  effect  of  the  interface  rough- 
ness, characterized  by  the  degree  of  lateral  disorder  77  given  in  (2). 
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Figure  3 shows  the  transmission  coefficient  calculated  for  different  values  of  ij.  Here 
Uy  = lOiiin,  (I,  = 0.3  nin,  = 20nin,  M = 50,  and  N = 38.  The  barrier  widths  are  2.1  nin 
for  both  the  emitter  and  the  collector,  their  heights  are  also  the  same,  0.3  eV,  and  the  well 
width  is  4.8  nm.  Thr('('  diffenuit  values  of  ij  were  studi('d,  namely  7/  = 0 (i)erfeet  DBH), 
Tj  = 0.3  nm  (largest  fluctuation  of  the  ord<-r  of  one  monolayer)  and  ?;  = 0.6  nni  (largest 
fluctuation  of  tin;  ord('r  of  two  monolayers).  As  a main  result,  it  can  be  seen  in  Fig.  3 
that  increasing  the  degree  of  lateral  disorder,  t],  results  in  a decrease  of  the  transmission 
probability  at  the  resonant  energy.  Notice  that  the  resf)iiant  p('ak  slightly  widens  due  to 
the  fluctuations  of  its  energy  for  each  reali/ation  of  the  disorder.  Besides,  an  additional 
effect  can  be  seen,  that  is,  as  the  degree  of  lateral  disorder  ;;  increases  the  conductance 
peak  ,shifts  to  smaller  energies.  The  lowering  of  the  energy  of  the  resonance  can  be 
utiderstood  a.ssuming  that  surface  roughness  makes  the  effective  width  of  the  ciuantum 
well  larger  than  its  nominal  value.  This  shift  woukl  make  both  the  threshold  voltage 
and  tin;  negative  differential  r('sist:ince  in  the  curnuit  voltage  characteristic  of  DBH  to 
appear  at  lower  bias.  But  an  effectiv('  wider  well  inipli(;s  a higher  current  too,  so  in  some 
statistical  sense  tlu'  current  for  the  disordered  DBH  will  l>e  higher. 
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Figure  3,  Transmission  coefficient  through  an  ordered  DBH  compared  with 
that  of  a disordc'red  DBH  for  two  different  values  of  the  degree  of  lateral  disorder 
]j.  The  disordered  results  comi)ris(‘  100  realizations  of  the  disorder. 


Regarding  tin'  effect  of  the  sizt'  of  the  islands,  ha\’e  observt'd  that  it  can  be 

neglected  unless  this  size  is  of  the  order  of  the  electron  wa\’ch'ngth,  that  is,  for  C K 
the  transmission  coefficient  does  not  depend  on  ((.  As  expected,  when  A,.  < C the  electron 
•starts  to  see  each  island  and  then  the  transmission  coefficient  decreases  as  increases,  as 
shown  in  Fig.  4.  For  energies  about  0.1  eV  this  transition  takes  place  at  sizes  ( ~ 10  nin  [3]. 
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Figure  4.  Transmission  coefficient,  through  imperfect  DBH  as  a function  of 
the  size  of  the  islands  for  ( ^ K.  (dotted  line)  and  C < (dashed  line).  Curves 
comprise  50  realizations  of  the  disorder.  Inset  shows  an  enlarged  view  of  the 
main  resonance  peak. 


CONCLUSIONS 

In  this  paper  we  have  presented  a method  to  study  electron  transmission  in  uninten- 
tionally disordered  heterostructures  with  rough  interfaces.  We  have  shown  that  the  main 
effects  of  the  lateral  disorder  are  to  decrease  the  transmission  trough  DBH  and  to  lower 
the  resonant  energy.  As  a consequence,  our  model  predicts  that  the  threshold  voltage  of 
resonance  tunneling  diodes  becomes  smaller  in  imperfect  quantum  devices. 
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